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Abstract. We study the interplay between the Kondo and Andreev-Josephson effects in a 
quantum dot coupled to one normal and two superconducting (SC) leads. In the large gap 
limit, the low-energy states of this system can be described exactly by a local Fermi liquid 
for the interacting Bogoliubov particles. The phase shift and the renormalized parameters for 
the Bogoliubov particles vary depending on the Josephson phase (j) between the two SC leads. 
We explore the precise features of a crossover that occurs between the Kondo singlet and local 
Cooper-pairing states as varies, using the numerical renormalization group approach. 



1. Introduction 

The Kondo effect in quantum dots (QD) connected to superconducting (SC) leads lias been 
studied intensively in these years, experimentally jUH] and theoretically [5Vll5j. The competition 
between the strong Coulomb repulsion and superconductivity gives rise to the quantum phase 
transition (QPT) between a nonmagnetic singlet and magnetic doublet ground states. It also 
causes the Andreev scattering at an interface between the QD and a normal lead, and affects 
substantially transport properties near the QPT, or the corresponding crossover region. 

In the present work, we consider a single QD coupled to one normal and two SC leads as 
illustrated in Fig. [TJ and study the interplay of the Kondo, Andreev, and Josephson effects, by 
using the Wilson numerical renormalization group (NRG) approach. We show that three phase 
variables, the phase shift 5 for the interacting Bogoliubov particles, the Bogoliubov angle Qb 
at the interface between the QD and normal lead, and the Josephson phase (j) between the two 
SC leads play a central role, and determine the low-temperature properties in the large SC gap 
limit. 
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Figure 1. Anderson impurity (•) coupled to one 
normal (A) and two superconducting leads {L,R): 
Fj, = npvl with p the density of states of the 
lead {u = L,R,N), and Vi, the tunneling matrix 

element. The SC gaps ^l/r = I^l/_rI e*^^/^ give 
rise to the Andreev scattering and Josephson effect. 
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Figure 2. The left panel shows the phase boundary between the non-magnetic singlet (upper 
side) and magnetic doulbet (lower side) ground states in the U vs Tg plane for several values of 
(p = — The middle and right panels show the Josephson current J and {d-^d^) as functions 
of (f> for r^. = 2. OA for several U. The critical current is given by Jc = eA/h for finite A. 



2. Model 

We start with a single impurity Anderson model of the form, 

H = ^dind-l) + ^ (rid -lf+ Yl (^i,-^- + '^-^-,-) + ^l-d + ^BCS > (1) 

u=N,L,R cr 

i^lcad= Yl 2^ ^fc4,fcaC^,fca, ^BCS= Y {^^ "^M 4,,-^ + ^■'^■) ■ (2) 

u=N,L,R k,a a=L,R k 

Here, = + U/2, U the Coulomb interaction, = X^^dJ-dcr) ^i^d dt (cli^^) the creation 
operator for an electron with energy (e^) and spin a in the quantum dot (lead). The QD and 
leads are coupled via u^, and V'jy^o- = c^^ ka/ V-^v with Mu the normalization factor one-particle 
states in the leads. We assume that = -Kpvl and the density of states p are independent 
of the frequency w. The complex superconducting gap A^^^ = |A^^^|e*^^/^ gives rise to the 
Josephson current between the two SC leads for finite <P = (^r — ^l- Since this Hamiltonian still 
has a number parameters to be explored, in this report we concentrate on the half-filled case 
Ed = —U/2, assuming a symmetric junction with T ^ = T (= r5'/2) and |A^| = |Ajj| (= A). 

3. Two SC terminal case without the normal lead r^r = 

We first of all consider a simpler F^r = case, where the QD is coupled only to the two SC leads, 
in order to see the competition between the Kondo and Josephson effects |10H14j . Figure [2] shows 
the NRG results for the ground state properties. We see that the region of the magnetic doublet 
state, which emerges for large U/A. and small F^/A, expands in the parameter space as the 
Josephson phase (f) increases. This indicates that (j) disturbs the Kondo screening. At the QPT, 
the Josephson current and the SC correlation {d^di^ at the impurity site vary discontinuously, 
and take small negative values in the magnetic doublet phase for finite A. 

4. Low-energy properties of the QD connected to one normal and two SC leads 

We next consider the three terminal case, F^v 7^ 0, with the additional normal lead. Specifically, 
we concentrate on the A — )• cxd limit, where A is much larger than the other energy scales 
A S> max(F5, Ftv, C^, |ed|), and the quasi-particle excitations in the continuum energy region 
above the SC gap are projected out. Nevertheless, the essential physics of the low-energy 
transport is still preserved in this limit, and a static pair potential A^ is induced into the 
impurity site owing to the SC proximity effects [TlfTO]. 

Ad = Vne^'n + TLe''^ = \Ad\ e''^ , |Arf| = F5 ^l-Tosin^ (0/2) , (3) 




Figure 3. NRG results for the phase shift 6 for the Bogoliubov particles (left), the conductance 
9ns between the QD and normal lead (middle), and the Josephson current J (right), are plotted 
vs (j) for several Tjv in the large gap limit. The other parameters are chosen such that U = l.Sr^, 
Tl = Tji (= Ts/2), and = —U/2. Here, Jc = eTg/h since A — oo. 



where To = AT ^/ {T ^ + F^)^. Carrying out the Bogoliubov transform with the angle 0^, 



cosGb = -^, sineB = ^, Ea = + |Arf|2, (4) 

the Hamiltonian can be mapped onto the asymmetric Anderson model for the Bogoliubov 
particles, the Green's function for which takes the form Giuj) = [w — Ea + iT^ — S(a;)]~^ [15j . 
The phase shift 5, defined by G(0) = — |G(0)|e*'^, determines the ground-state properties 

[n^) - 1 = (25 /tt - 1) cos , {d^d^) = ^ {26 /-k - 1) e'^d sin 9^ . (5) 

The conductance g^qg between the QD and normal lead, and the Josephson current J between 
the two SC leads can be determined by the three phase variables 6, Qb and (f), 



4e2 . 2 „ . 2 OA T To \26/'K-l\ sinG^ sin, 

— — sm G^ sm 26, J = — , 

/i h 2Vl-rosin^ {(t)/2) 



Furthermore, the renormalization factor Z = [1 — dTi/doj] ^, the position of the Andreev level 
Ea = ^Z[Ea + S(0)], and the Wilson ratio R for the Bogoliubov particles can also be deduced 
from the local Fermi-liquid behavior, using the NRG. 

The results are plotted vs (j) in Figs. [3] and H] for several values of F7v/r5. The Coulomb 
interaction is chosen to be f/ = l.SF^, and for this value of U the QPT occurs in the F^r — t- 
limit at (/) ~ 0.467r, where Ea = U/2. Note that To = 1 and = in the present case. The 
gapless excitations near the normal lead change the sharp QPT into a continuous crossover 
between the two different singlet states. We see in Fig. [3] that for cf) < 0.467r the ground state 
is a local Cooper pairing consisting of a linear combination of the empty and doubly occupied 
impurity states with a small 6. On the other side, for > 0.467r, the ground state is a Kondo 
singlet with 6 ~ 7r/2. At the transient region, the conductance g^^g has a peak, and the 
Josephson current J decreases rapidly to zero. Note that a weak current, which flows in the 
opposite direction as seen in Fig. [2] for the magnetic doublet state with finite A, is absent in the 
large gap limit. These features that are caused by the crossover are smeared as F^r increases. 

Figure m shows the results for the renormalized parameters. We see for small F^r (= O.OSF^) 
that the parameters are strongly renormalized in the Kondo-singlet region for (j) > 0.467r, where 
Z <^ 1.0, R ~ 2.0, and the couple of the Andreev peaks Ea stay close to the Fermi level a; = 0. 
In contrast, in the local Cooper-pairing region for (p < 0.467r, the parameters are almost not 
renormalized Z ~ 1.0, i? ~ 1.1, and the Andreev peaks situate away from the Fermi level. When 
the coupling Tjy between the QD and normal lead is large, these two singlet states cannot be 
distinguished clearly. 



Figure 4. NRG results for the level position of renormalized Andreev resonance Ea (left), the 
renormalization factor Z (middle) and R — 1 (right), where R is the Wilson ratio, are plotted 
vs (j) in the large gap limit for U = l.5Ts and Tl = Tji (= Ts/2) for several Fat. 

5. Summary 

We have studied the crossover between the Kondo singlet and the local Copper pairing, which 
consists of the linear combination of the empty and doubly occupied impurity states, in the 
quantum dot coupled to the one normal and two SC leads. In this three terminal geometry 
the normalized parameters, which characterize the Fermi-liquid behavior of the Bogoliubov 
particles, vary depending on the Josephson phase (f>. Therefore, the crossover can occur at finite 
Josephson phase (j)Q. We calculated the phase shift 6 and the renormalized parameters with the 
NRG approach, and observed numerically that the Andreev conductance between the dot and 
normal lead has a peak near the crossover region at ~ (pc- The Josephson current between the 
the two SC leads decreases rapidly in the strongly renormalized Kondo-singlet ground state for 
(p > (j)Q, where the Andreev states Ea stay close to the Fermi level. In the local Cooper-paring 
ground state at the opposite side 4> ^ 4>(j-, the local Fermi- liquid parameters are not renormalized 
so much. 
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